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ABSTRACT 

Recent observations and detections of interstellar objects (ISOs) passing through the solar system have sparked a wave of interest 
into these objects. Although rare, these ISOs can be captured into bound orbits around the Sun. In this study, we investigate the 
novel idea of capture of ISOs into near-Earth orbits and find that a steady population of ISOs exists among the current population 
of Near Earth Objects (NEOs). Using numerical simulations, we find that the capture of ISOs into near-Earth orbits is dominated 
by Jupiter which is 10*x more efficient in capturing ISOs. Captures are more likely to occur for objects with high eccentricities 
and low inclinations. We also investigate the stability of captured ISOs and find that they are generally unstable and survive 
shorter than known NEOs with a half-life time of ~ 0.05 Myr and are ejected from the solar system due to interactions with other 
planets or the Sun. Our results have important implications for understanding the population of interstellar objects in the solar 
system and possible future detection. We find that about 1 — 2 50 — 70 m sized captured ISOs among NEOs would be detectable 
by LSST over its lifetime. By detecting and studying captured interstellar objects, we can learn about the properties and origins 


of such objects, and the formation and evolution of exoplanetary systems and even our solar system. 
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1 INTRODUCTION 


The discovery and identification of interstellar objects (ISOs) 
1I” Oumuamua (e.g., Meech et al. 2017; Micheli et al. 2018) and 
2I/Borisov (e.g., Jewitt & Luu 2019) recently has led to a flurry of 
interest in the origins and dynamics of these objects. Siraj & Loeb 
(2022b,c) also reported the discovery of two apparent interstellar me- 
teoroids CNEOS 2014-01-08 (IM1) and CNEOS 2017-03-09 (IM2) 
which produced fireballs that were detected by satellites. IM1 was 
later confirmed to be of interstellar origin by U.S. Department of 
Defense (Siraj & Loeb 2022b). Several hypotheses have been put 
forward to explain the dynamical origins of such objects. For exam- 
ple, Portegies Zwart et al. (2018); Portegies Zwart (2021) and Torres 
et al. (2019) considered the possibility of interstellar interlopers as 
objects that had been ejected from the debris disk due to encounters 
with other stars in the birth cluster of the star. Moro-Martin (2018, 
2019) also considered the possibility that such objects originated in 
young protoplanetary disks or were a part of an exocometary cloud. 
Hands et al. (2019) explored the exchange of these objects between 
different stellar systems in an open cluster environment. Cuk (2018) 
found that the origins of such objects were better explained as a tidal 
disruption fragment from a binary star system. It is, therefore, imper- 
ative to study ISOs as they can open up avenues for understanding 
origin and dynamics of objects beyond the solar system (e.g., Selig- 
man & Laughlin 2020; Hoang & Loeb 2020; Siraj & Loeb 2022a; 
Bergner & Seligman 2023; Hoang & Loeb 2023). 
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Although ISOs are rare, their discovery suggests that they could 
be captured by planetary systems. The capture of interstellar objects 
by the solar system has been studied for a long time (e.g., Valtonen & 
Innanen 1982; Torbett 1986). Most studies have explored the cross- 
section of capture of the Sun-Jupiter system and captures as a result 
of close encounters with the giant planets (e.g., Siraj & Loeb 2019; 
Hands & Dehnen 2020; Napier et al. 2021a; Dehnen & Hands 2022; 
Dehnen et al. 2022). Jupiter is expected to dominate the capture 
cross section and the rate of capture in the present day solar system is 
expected to be low due assuming that the hyperbolic excess velocity 
of the ISOs is similar to the field star velocity distribution (e.g., 
Napier et al. 2021a; Dehnen et al. 2022). However, the capture of 
ISOs would be easier when the Sun was still residing in its birth 
cluster given the low velocity dispersion of field stars in the birth 
cluster. 

The question of whether any captured ISOs reside in the solar 
system currently is of major interest and is a matter of investiga- 
tion. Morbidelli et al. (2020) report that no known planetesimals of 
interstellar origin currently reside in the solar system. Napier et al. 
(2021b) simulated a suite of 270k synthetic captured objects and 
found that most captured objects were not able to raise their perihe- 
lion efficiently leading to repeated scatterings with the giant planets 
which make them unbound over time. Only 3 out of 270k objects 
simulated survived for a duration of 1 Gyr. Therefore, to constrain 
any useful information about any possible population of interstellar 
objects in the current solar system, more simulations are required. 

Due to the larger number density of smaller interstellar objects, 
we would except smaller captured ISOs to outnumber any larger 
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ISOs, should they be present. The smaller objects would be harder 
to detect unless they were close to Earth. This compels us to ask an 
interesting question: does a population of ISOs lie amongst the Near 
Earth Objects (NEOs) today? Although the capture rate of ISOs is 
expected to be dominated by the giant planets, especially Jupiter, 
Earth can play a role as well in capturing ISOs into NEO orbits. An 
interesting but unexplored avenue of exploration is related the capture 
of ISOs by the Earth-Moon system. We expect the number of smaller- 
sized ISOs to be much larger than the number of ?Oumuamua and 
Borisov sized ISOs (e.g., Siraj & Loeb 2022b). This would ensure 
that despite the smaller capture cross-section, the Earth-Moon system 
would be able to capture some of these smaller sized objects. The 
orbital parameters of such a population would need to be identified 
and compared to that of present day NEOs. If a population of such 
captured objects exist among present day NEOs, they could be more 
easily detectable and can be visited by space probes which would 
give us invaluable information about the formation and evolution of 
exoplanetary systems and even our own solar system. This compels 
us to study whether interstellar NEOs exist and therefore motivates 
this work. 

In this study, we explore the dynamics and capture of objects by 
the Earth-Moon and Jupiter systematically using a large suite of N- 
body scattering simulations. The total number of scattering events 
simulated is ~ 10!! which can be considered state-of-the-art. This 
is required given the small capture cross-section of the Earth-Moon 
system. We also explore the role of the Moon on capture of ISOs by 
the Earth. In addition, we compare the dynamics of captured objects 
by the Earth-Moon system with that of a giant planet like Jupiter and 
explore whether close encounters with Jupiter can lead to production 
of any NEOs. After performing long-term simulations, we discuss 
the possibility of finding ISOs among a population of NEOs today. 

We begin by introducing a new hybrid integration scheme that 
enables us to perform a large number of simulations in Section 2. 
This is followed by a discussion of our models and initial conditions 
in Section 3, and results in Section 4. We, then, discuss our results in 
the context of previous works and compare our captured population 
to the population of small bodies in the Solar System at present day 
in Section 5. We conclude in Section 6 and talk about some future 
work that needs to be done to further this research. 


2 METHODS 


In order to derive robust statistical results, a large number of scat- 
tering experiments need to be performed. A full numerical solution 
involves numerically integrating the orbit of the planet (Earth/Jupiter) 
and any associated Moons taken into consideration and the orbit of 
the incoming ISO. Since the trajectory of the incoming ISO is far 
from the orbit of the planet at most times, a full numerical solution 
becomes quite inefficient since the orbits of the planets are approxi- 
mately Keplerian. This presents a computational bottleneck. In order 
to alleviate this issue, we have derived a new hybrid method that 
is much faster than full numerical solution while producing similar 
results as that of the latter. 

Since the scatterings occur over a time-scale that is much shorter 
than the secular timescale of the system and the incoming test particle 
has zero mass, we can model the orbit of Earth-Moon/Jupiter around 
the Sun analytically as a pure Keplerian orbit. In order to do it, we 
split the system hierarchically into two subsystems. For Earth-Moon, 
the movement of the Earth-Moon barycenter around the sun is one 
subsystem whereas the movement of the Moon around the Earth 
is another subsystem. For Jupiter, we don’t consider any associated 
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moons and thus we only evolve its orbit. For the timescales in question 
for the scattering experiments, we compare the motion of the system 
using our hierarchically split Keplerian method to those from full 
numerical results and find similar results. The positions and velocities 
of the particles in the system are generally agree within 1078 — 19-10 
of each other. We use a Kepler solver that uses the universal variable 
formulation (Wisdom & Hernandez 2015). 

To integrate the motion of the test particle under the influence of 
the Sun-Earth-Moon system (or any other system for that matter), 
we use a time-symmetrized fourth-order Hermite scheme(Kokubo 
et al. 1998). For more information, we refer the interested reader to 
Kokubo et al. (1998). 

We perform long-term integrations for a subset of captured objects. 
For these simulations, we include every planet except Mercury to 
understand the lifetime of captured ISOs in NEO orbits. We again use 
the Kepler solver to integrate the orbits of the planets without taking 
into account any interactions between the planets. The Keplerian 
portion of the integration is done in Jacobi coordinates (e.g., Murray 
& Dermott 2000; Rein & Tamayo 2015). Essentially our integrator 
then reduces to a Wisdom-Holman (Wisdom & Holman 1991) like 
scheme without the interaction terms between planets taken into 
account. Although we expect the interaction terms to contribute to 
the secular evolution of the planets, we expect the effect of the secular 
evolution of the planetary orbits on the orbit of the ISO to be minimal. 
We stress that the usage of the Hermite integrator along with the 
Keplerian scheme allows us to simulate close encounters which are 
necessary for the long term simulations. 

We time-symmetrize the integrator by using an iterative P(EC 3 
scheme (Kokubo et al. 1998) along with symmetrized timesteps as 
used in Pelupessy et al. (2012). We use a globally adaptive timestep 
where the timesteps are determined only by the scattered particle. 
Time steps are taken as the minimum of the following timesteps: 


Fij 


J 
Tfreefall = 17 ae (1) 
1J 
Thyby = n7 (2) 
yoy =N 


where r;j is the relative distance, v;; is the relative velocity, and 
aij is the relative acceleration between two particles. 7 is tunable 
parameter that controls the accuracy of integration. We set it to 1073 
in our simulations. 

The Hermite-Kepler hybrid integration scheme proves to be ex- 
tremely efficient at simulating the scattering of ISOs and their long 
term evolution. We compare our results to those obtained using the 
IAS15 integrator (Rein & Spiegel 2015) for a set of 107 scattering 
events. We notice that the our method provides statistically simi- 
lar results as that from IAS15 (see appendix). The experiments are 
~ 10 — 50x faster when using this method compared to IAS15 en- 
abling a larger number of scattering experiments to be performed. 
Our hybrid integrator also offers the ability to handle close encoun- 
ters unlike a traditional Wisdom Holman scheme, allowing us to 
examine the long term evolution of captured ISOs. 


3 MODELING THE SCATTERING OF OBJECTS 


To model an incoming interstellar object on an initially unbound or- 
bit, we follow the steps outlined in Napier et al. (2021a) and Dehnen 
& Hands (2022) with some differences. We consider scattering events 
with Earth-Moon and Jupiter. The process of generating initial con- 
ditions is similar in both cases and we present an overview of our 
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Figure 1. Visualization of trajectories of incoming particles getting scattered by the Sun-Earth-Moon system. Since there are two binary systems present, we 
can calculate the cross section of the whole system and that of the Earth-Moon system. Left: A particle experiencing a close encounter with the Sun. Right: A 
particle experiencing a close encounter with the Earth-Moon system. Image credit: NASA. 


methods using the example of Earth-Moon before moving on to 
Jupiter. 

In the first scenario, our system only consists of the Sun, and the 
Earth-Moon system. This is done in order to study the capture cross 
section and the influence of only the Earth-Moon system without 
that of the giant planets. We first explore the cross-section of the 
entire Sun-Earth system before moving on to simulate trajectories 
that undergo close encounters with the Earth-Moon binary. The initial 
conditions for the incoming particles are generated in different ways 
in the two scenarios as is detailed below. 

Because of the time dependent potential of the Earth-Moon system, 
ISOs that pass within the orbit of the Earth around the Sun can be 
scattered and captured. This determines the cross-section of capture 
for the Sun-Earth system. In order to calculate this cross-section, we 
follow steps as outlined in Napier et al. (2021a) to generate the initial 
conditions. 

As in Napier et al. (202 1a), the incoming particle is sampled from 
a sphere with a radius of 10° au from the barycenter of the Sun-Earth- 
Moon system. Since we are interested exploring the cross-section as 
a function of the hyperbolic excess velocity væ, we simulate multiple 
scattering events for a particular vo. Once vo is chosen, the speed 
of the particle at the current radius (v,) can be derived using the 
conservation of energy as 


2G 
ra 


3 
dary (3) 


vr = 
where u is the barycenteric mass of the Sun-Earth-Moon system and 
dhary İs the distance of the particle from the barycenter of the system. 
The velocity vector for the particle is then found by multiplying 
the velocity from equation 3 with the direction directly towards the 
barycenter of the system. 

Given a choice of væ, the maximum impact parameter (bmax) of 
the object can be derived. The maximum impact parameter is given 
as 


bmax = y dha + 2¢max|a| (4) 


where gmax is the maximum perihelion distance sampled and a is the 
semi-major axis of the incoming particle. We set qmax = 1.3 au in 
our simulations. This is mainly motivated by the definition of NEOs 
in literature where objects with perihelion < 1.3 au are considered to 
be NEOs (e.g., Morbidelli et al. 2002, and references therein). The 
impact parameter of the object is sampled uniformly between 0 and 
bmax and the object is placed on tangent plane to the sphere at the 
point in arandom direction. Some previous studies have sampled the 
impact parameter assuming a probability distribution œ b°? (Siraj & 
Loeb 2019; Hands & Dehnen 2020; Dehnen & Hands 2022). In such 
a scenario, more encounters take place closer to the planet rather than 
far from it. We ran a set of simulations assuming the distribution of 
impact parameters was œ b? and found no major differences in the 
computed cross-section. 

To save on computation time, we then evolve the motion of the 
particle and the barycenter of the system as a pure Keplerian orbit 
until the particle has reached a distance of 20 au from the barycenter. 
This value is similar to the initial starting distance used in Dehnen & 
Hands (2022) but lower than that in Napier et al. (2021a). We verified 
that the lower value used in our simulations do not affect the results. 

Once the particle has reached a distance of 20 au, we evolve the 
particle under the influence of the Sun and the Earth-Moon system. 
We put in the Earth-Moon system around the Sun. The positions 
of the Earth and Moon are generated with random orbital phases 
with the orbital phase being drawn uniformly between 0 — 27. This 
is done to prevent any spurious results that may arise out of any 
specific configuration of orbital phases. The system with the test 
particle is then evolved until the particle has been scattered off by the 
system. In case the specific energy of the particle drops below zero, 
we stop the simulation and store the positions and velocities of all 
particles in the simulation. In case of collision with any of the other 
objects in the system, we also record the positions and velocities of all 
particles in the system. If the outgoing test particle is unbound after 
it reaches a distance > 20 au after scattering, we stop the simulation 
and conclude that the particle could not be captured. For all of the 
captured particles, the orbit is evolved until the particle has reached 
a distance > 20 au and then the orbital parameters are calculated. 
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To generate trajectories that undergo close encounters with the 
Earth-Moon, we follow the steps above and then select orbits that 
undergo close encounters with the Earth-Moon system, similar to 
Dehnen & Hands (2022) with some differences. We define a close en- 
counter when the minimum distance between the incoming ISO and 
the Earth is lower than than the Hill radius of Earth (0.01 au),similar 
to within L1 and L2. 

To select the subset of systems wherein the incoming ISO un- 
dergoes close encounter with the Earth-Moon, we chose only those 
orbits where the hyperbola of the orbit of the incoming ISO inter- 
sects with the orbit of the Earth-Moon system and lies within the Hill 
radius of Earth. This is a non-trivial problem as it involves finding 
the Minimum Orbit Intersection Distance (MOID) and selecting a 
subset where the MOID is less than 0.01 au. To accomplish this, we 
use a fast, geometric method to find the MOID using the technique of 
Wizniowski & Rickman (2013). The calculation is accurate to 10714 
au and any missed MOIDs do not affect our results. 

Once an orbit has been found, we evolve the particle on this hy- 
perbola, and place the Earth-Moon system at a random orbital phase 
such that the incoming ISO is at the MOID distance with respect to 
Earth at t = 2 yr. The orbit of the ISO is simulated for a total for 
t = 5 yr after which the specific energy of the ISO is evaluated to 
check whether it was captured or not. Collisions with the Earth-Moon 
system are also tracked. An additional set of simulations without the 
Moon is performed using the same initial conditions. This is done in 
order to study the effect of the Moon on captures by close encounters 
with the Earth-Moon system. 

We also perform a set of 10? simulations with the Sun-Jupiter 
system in order to understand how the orbital properties of the cap- 
tured object differs when captured by the Earth-Moon system versus 
that by Jupiter and to study whether Jupiter can efficiently capture 
interstellar NEOs. The process of generating initial conditions for 
ISOs scattered by Jupiter is similar to that mentioned above. For the 
Sun-Jupiter system, we perform a full integration once the particle 
has reached a distance < 100 au. qmax is set to 7.3 au for these 
Sun-Jupiter simulations. 

To study the lifetime of the captured objects, we perform long- 
term simulations upto 10 Myr using the hybrid scheme described in 
the previous section. Once the scattering simulations have ended, we 
place the captured object and the planetary system at their respective 
positions in the solar system. Our solar system consists of all planets 
except Mercury whose effect is expected to be minimal. Since the 
long-term simulations are computationally expensive even with the 
hybrid integration scheme, we restrict ourselves to performing inte- 
grations of ~ 20k captured objects with a < 20 au that are randomly 
selected. 


4 RESULTS 
4.1 Capture cross section of the Sun-Earth-Moon system 


We performed a total of 5 x 10!9 simulations to calculate the cross- 
section of capture of the Sun-Earth system. More simulations were 
performed for simulations with higher voo because of the smaller 
capture cross section. To calculate the cross-section of capture, we 
used the following formula (Dehnen & Hands 2022): 


N captured = y N, captured 


N sampled 


o= nb? ax (5) 


To calculate Ncapturea, We only consider particles with semi-major 
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axes a < 10° au. Particles with larger values of a can be stripped 
by nearby stars and become unbound from the system. We also 
investigated using a different formula for calculating o using equation 
(29) from Napier et al. (2021a) but found minimal differences. 

Looking at Figure 2 where we plot the capture cross section (a) 
as a function of væ for both Earth-Moon and Jupiter, we find that o 
behaves as vo at lower values of vœ and as v at higher values. 
This is qualitatively similar to the results from Napier et al. (2021a). 
We fit a function of the form 


90 


(s) (sy 
Vo Vo 

to the measured datapoints and find that og = 9.825 x 104 + 3.55 x 
10%au? and vp = 0.0914 + 0.0011km/s. The latter of these values 
represents the velocity at which the capture cross section transitions 
from a v% form to a vo form. We calculate the cross-section of 
Jupiter in a similar fashion by performing 10° scattering simulations 
and compare it to that of Earth-Moon system. Although they follow 
the same functional form, Jupiter transitions from a v% form to a vio 
form at higher væ. Indeed we find that vo jup = 0.554 + 0.020km/s. 


We find that 7 jup = 6.09 x 104 + 4.47 x 103 au”. The cross-section 
of capture for the Sun-Earth system is 1% of that of Jupiter at smaller 
Voo and 0.01% that of Jupiter at larger væ. Although our fitted values 
are slightly different from that found by Napier et al. (2021a), we 
note that our fitted cross-section only differs by < 10% compared 
to theirs. This difference can be attributed to the presence of other 
giant planets like Saturn while calculating the capture cross section in 
Napier et al. (2021a). However, the similarity in the results pertaining 
to the capture cross section of Jupiter validates the robustness of our 
computational methods. 

Objects that are captured can have varying periastron distances 
with respect to the planet. We define weak encounters as those where 
the minimum distance of the incoming object to the planet is greater 
than the Hill radius of the planet and close encounters as encounters 
where the minimum distance is within the Hill radius of the planet. 
Both weak and close encounters can result in the capture of ISOs. 
To understand which type dominates as a function of the hyperbolic 
excess velocity, we plot the distribution of particles as a function of 
the minimum distance to Earth dg min and the semi-major axis a 
of the captured objects in Figure 3 We find that at lower vo. weak 
encounters dominate the capture. Weak encounters, in fact, account 
for ~ 99.7% of the captures at low hyperbolic excess velocities. The 
ISOs captured due to weak encounters tend to have larger semi- 
major axes and smaller inclinations. The median semi-major axis of 
objects captured due to weak encounters at Voo = 0.1 km/s is 6 X 104 
au, whereas that for close encounters is 2.2 x 10° au. The median 
inclination for weak encounters is 59° while that for close encounters 
is 98°. At voo > 1 km/s, close encounters become the dominant 
mechanism to capture ISOs. Close encounters permit captures of 
objects with larger inclinations and at smaller væ even objects in 
retrograde orbits can be captured. 

Next, we analyze the orbital parameters of the captured objects 
in Figure 3. As væ increases, looking at Figure 3, we find that the 
distribution of the semi-major axis of captured objects moves towards 
lower values. The median of the distribution at væ = 0.1km/s is 
6 x 104 AU whereas that for Voo = 2.0km/s is about 5.5 x 102 AU. 
This leads to a fatter tail in the eccentricity distribution as well. At 
higher væ most captures take place on orbits that are prograde with 
respect to the orbit of the Earth around the Sun. Orbits on lower 
inclinations allow for longer interaction time between the incoming 
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Figure 2. Capture cross section (o`) presented as a function of the hyperbolic excess velocty (vo). We have compared the cross section of the Earth-Moon to 
that of Jupiter. The dots represent the actual data points calculated from the simulation and the lines represent the fitted cross section. We note that the cross 
section of Earth-Moon follows the same qualitative form as that of Jupiter. However, the cross-section of the Earth-Moon system is 1% of that of Jupiter at lower 


velocities whereas it is 0.01% that of Jupiter at higher velocities. 


ISO and the Earth-Moon system leading to easier captures. There 
is a pronounced peak at i = 90° but this is caused due to the initial 
distribution of inclinations in our system. 


4.2 Close Encounters with Earth-Moon System 


The fraction of orbits that undergo close encounters with the Earth- 
Moon is minuscule in the simulations presented in 4.1. As we saw 
in the previous section, close encounters are the dominant source 
in production of captured objects at larger vo. and objects that have 
smaller a. Since we are interested in objects that have been captured 
into the solar system (a < 30 au), we expect close encounters with 
the Earth-Moon to be the dominant mechanism rather than weak 
encounters. Here we present a closer examination of the captures 
resulting from close encounters and the cross-section of capture into 
smaller semi-major axis. 

Since close encounters occur less frequently, the method of pre- 
vious section proves to be extremely computationally expensive to 
obtain better statistics on objects captured due to such encounters. 
Therefore, we perform a different set of simulations where the orbit 
of the incoming ISO falls within the Hill radius of the Earth using 
the method described in section 3. We perform a set of 4 x 10!9 sim- 
ulations to calculate the cross-section of capture of the Earth-Moon 


system, the orbital parameters of the captured objects, and the effect 
of the Moon on the capture cross-section of the Earth. 
Only a fraction fe, where fe is given as 


r2 
Hill 
fe = (7) 
Imax 
of the total scattering events result in a close encounter. Therefore, to 
calculate the capture cross-section resulting from close encounters, 
we modify Equation (5) slightly to obtain 


N, captured + 4y N, captured 


(8) 
N sampled 


Oc = fenbhax 
where Nsamplea is now the total number of close encounters sampled. 
To test the validity of this method, we compare the cross section of 
capture at larger Vo from the previous section to the cross section 
obtained for objects captured with a < 10° au using this method. 
Since the captures at larger vœ are dominated by close encounters, 
we should expect both methods to give us consistent values for the 
capture cross-section at larger Vo. We find that both methods produce 
results that are within 1% of each other indicating this method to 
calculate the capture cross-section of close encounters is valid. 

We present the cross section of capture due to close encounters 
(Tc) calculated using equation 8 in Figure 4 in order to understand 
how the cross section changes as a result of changing the threshold 
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Figure 3. Distribution of the orbital parameters of the objects captured by the Sun-Earth system. Top: Distribution of the minimum distance to Earth dg min 
and the captured semi-major axis a. The brighter colors denote larger number of captued objects whereas darker colors indicate fewer. At vo = 0.1 km/s, the 
dominant source of captures are wide encounters. As Væ increases, more captures occur because of close encounters with the Earth-Moon system. In fact, all 
captures for Vo > 1km/s occur due to close encounter with the Earth-Moon system. Bottom: Semi-major axis (a), eccentricity(e), and inclination (i) distribution 
of captured objects. As væ increases, the median of the semi-major axis distribution decreases. We also find that prograde orbits are preferred for captures. 


semi-major axis. We calculate the cross section for three different 
values of the threshold semi-major axis a;: 10 au, 20 au, and 50 au. 
We only select particles with semi-major axes below the threshold 
semi-major axis to calculate Ncapturea- We find that the cross-section 
falls off as v32 until 10 km/s and then it drops rapidly. 

Fitting a functional form to the capture cross section of close 
encounters dc is more complicated than that presented in the previous 
section for all encounters. In order to do so, we assume that the cross 
section drops off as v5® after vo = 15 km/s. Note that while the 
assumption might not exactly hold, it makes little practical difference 
since the capture cross section is dominated by smaller vo where the 
functional form is well known. As in the previous section, we fit a 
functional form similar to equation 6. However, oo and vg are now 
functions of the threshold semi-major axis a; and written as 09 (ar) 
and vg(a;). We find a weak dependence of vg as a function of a; 
where vo(ar) decreases as a function of a+. However, the same is not 
true for o9(a;). Comparing o(a;) as a function of a+, we find that 
it can be best represented as 


colar) = Tora? (9) 


. We find the best fit values to be oo; ~ —1.1x 107!9+1.5x 107 au? 
and £ = 3.07 + 0.10. Thus, we find that ce (ar) œ a?-07+0.10, How- 
ever, we caution the reader this propotionality does not include the 
contribution from vg(a;z). 
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For væ > 10 km/s, the capture cross-section drops rapidly. At væ = 
10 km/s, the cross-section of capture is ~ 107°au?. The geometric 
cross-section for collision with the earth is Geo) ~ 5 x 10-2au2. 
Thus, for væ = 10 km/s or higher, collisions with the Earth dominate 
over captures. 

The distribution of orbital parameters is presented in Figure 5 and 
presents an interesting story. We notice qualitatively similar trends 
as that from Figure 3 but we notice that the distribution of semi- 
major axes of captured objects skews towards smaller values than 
those from wide captures. As væ increases, the distribution skews 
towards smaller values of a and in fact for væ = 8 km/s, the median 
value is 50 au. This is on account of closer approaches taking place 
with the Earth as is evident from the distribution of distances of 
closest approach to Earth (rp). The distribution of inclination of the 
captured objects shows that with larger values of vo, the inclination 
skews towards i = 0°. There is a smaller fraction of captured objects 
in retrograde orbits for væ = 0.2, 2 km/s but this drops rapidly as væ 
increases and we find no retrograde captures for vo > 8 km/s. 

An interesting examination pertains to the effect of the Moon on 
captures by the Earth. We run a set of simulations with the same 
initial orbits but with the Earth without the Moon to calculate the 
cross-section. We select objects with semi-major axes below a certain 
threshold and vary this threshold parameter. For all values of the 
threshold parameter > 10 au, we find little differences between the 
capture cross-section of the models with the Moon versus those 
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Figure 4. The capture cross section of close encounters (oc) with Earth-Moon system as a function of the hyperbolic excess velocity (væ) for three different 
sets of maximum semi-major axes(a). Since close encounters are rare, the capture cross section is much lower than the full capture cross section calculated in 


Figure 2. However, the capture cross-section follows the same qualitative functional profile as that in Figure 2. 


without. However when we go to a threshold semi-major axis of 5 
au, we find certain differences at higher væ. We find that the Moon 
actually has very little effect on the capture-section of the Earth until 
Væ ~ 12km/s. However, past that value, the capture cross section of 
the Earth-Moon system increases rapidly compared to the Earth-only 
system. At Væ = 13 km/s, the ratio of the capture cross section with 
the moon to the capture cross section without the moon is 1.3. At 
Voo > 13 km/s models with the Moon can capture objects with a < 5 
but models without the moon are unable to. At larger voo, the time 
dependent potential of the Moon around the Earth plays a larger role 
in the capture of ISOs and thus leads to more captures. This effect, 
however, is practically negligible as the captures are dominated by 
captures at lower vo. where the differences are negligible. 


4.3 Captures by Jupiter and Comparison to Known Small 
Bodies 


Since Jupiter is expected to dominate the capture cross-section, we 
examine the effectiveness of Jupiter in capturing objects into near 
Earth orbits. For the sake of clarity, we define an object to be in 
a near Earth orbit if its perihelion distance q is < 1.3 au (e.g., 
Morbidelli et al. 2002). We run a set of 10° simulations with the 
same væ as that used while simulating the Earth-Moon system and 
we set max = 7.3 au. We generate close encounters in a similar 
fashion as that mentioned before. We are primarily interested in 
close encounters with Jupiter. We calculated the efficiency of weak 
encounters in generating captured objects in the solar system and 
found that weak encounters become subdominant for large vo. and 
a < 100 au. To capture objects on orbits with a < 50 au, a close 
encounter with Jupiter is required. In this study, we consider all 


trajectories that pass within the Hill radius of Jupiter (0.35 au) as a 
close encounter. 


To understand Jupiter’s efficacy at injecting captured objects into 
NEO orbits, we consider objects with a < 10 au and q < 1.3 au 
and calculate the cross-section in Figure 6. Examining the figure, we 
find that the capture cross-section of close encounters with Jupiter 
qualitatively follows the same form as that of close encounters with 
the Earth-Moon. However, Jupiter, owing to its mass, is 104x more 
efficient than Earth-Moon in generating NEOs from ISOs. This is 
not surprising since the capture cross section for planetary close 
encounters (Frc) is « mi, where mp is the mass of the planet (Napier 
et al. 2021a, equation 22). Jupiter is ~ 300x more massive than Earth 
so we expect Ler. = 300° ~ 104. Performing a similar analysis as 


c,earth 
that in previous sections, we find that for Jupiter 7 (az) œ az 0. 10 


and p,e ~ (1.35 + 0.40) x 10~6au?. For Jupiter, we also find a 
dependence of vo(at) on az. It is given as volat) = ((15.44 + 
1.10)km/s)a70-36+0-02, 


Notably, objects captured by Jupiter have a different distribution 
of orbital parameters than those captured by Earth. In Figure 7, we 
calculate the distribution of captured objects by Earth-Moon and 
Jupiter at vo = 8 km/s and compare those to present day known 
NEOs. We choose a larger væ for presentation since the hyperbolic 
excess velocity of any incoming ISOs in the present day solar system 
is expected to be large. The median semi-major axes of captured 
objects with a < 100 au by the Earth-Moon at væ = 8 km/s is 32 au, 
whereas that for objects captured by Jupiter is 16.5 au. Jupiter is also 
more efficient at capturing objects on retrograde orbits. In fact, the 
Earth-Moon system is unable to capture any objects on retrograde 
orbits at higher væ. Jupiter is also able to capture certain objects 
on lower eccentricities than the Earth-Moon system. The minimum 
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Figure 5. Orbital properties of captured objects undergoing close encounters with the Earth-Moon system and their relative fractions. Top-left: Distribution 
of semi-major axis (a) Top-right: Distribution of pericenter distance from Earth (rp) Bottom-left: Distribution of eccentricity (e) after capture. Bottom-right: 
Distribution of inclination (i) after capture. We notice that as the hyperbolic excess velocity væ increases, the median semi-major axis a and eccentricity e of 


the captured objects decrease. Most objects are preferentially captured on prograde orbits and objects with higher velocities need closer encounters with Earth 
to get captured. The distributions are qualitatively similar to those in Figure 3. We find that no captures on retrograde orbits are possible for large væ. 


eccentricity of an object captured by Earth-Moon is 0.76, whereas 
that by Jupiter is 0.66. 


How do the orbital parameters of objects captured compare to 
the known population of NEOs today? To examine that, we plot the 
orbital properties of the known NEOs and the orbital properties of 
objects captured by Earth-Moon and Jupiter in Figure 7. Just for 
clarity, we have only considered a specific set of objects captured 
with a < 100 au at vœ = 8 km/s. We find that NEOs have a median 
semi-major axis of 1.7 au and are heavily skewed towards prograde 
and circular orbits. Most NEOs originate within the orbit of Jupiter. 
However, there is a smaller albeit significant fraction of NEOs orig- 
inating from a > 10 au. Our simulations cannot generate any NEOs 
with a < 2 au. The object with the minimal semi-major axis cap- 
tured by Earth is 4.5 au whereas that for Jupiter is 2.1 au. Objects 
cannot be captured on semi-major axes smaller than these values as 
that requires an encounter with the planet with periapsis smaller than 
the radius of the planet, resulting in a collision rather than capture. 
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Comparing the distributions, we find that, should an ISO were to 
be captured into an NEO orbit, it would more likely end up with 
a > 10 au. Incidentally, this is where the Centaurs exist. ISOs hiding 
amongst the Centaurs have been examined by Siraj & Loeb (2019) 
but no known Centaurs are considered to have an interstellar origin. 
However, our study suggests a closer examination may be merited. 
We caution the reader that this is not a full representation of present- 
day captured ISOs by Earth and Jupiter (should they exist). ISOs 
captured in the past will evolve their orbital parameters over time due 
to interactions with other planets in the Solar System. ISOs captured 
in the present day would also arise from a different velocity distribu- 
tion which can affect the distribution of orbital parameters. We defer 
the presentation of evolution of orbital parameters to future work. 
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Figure 6. Capture cross section of close encounters (oc) of objects with a < 10 au and q < 1.3 au as a function of the hyperbolic excess velocity (voo). We 
note that rate of capture falls off as vg? at smaller vo and drops rapidly at larger væ. Jupiter is 104 times as efficient as Earth in injecting objects into the NEO 
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Figure 7. Distribution of known small bodies in the solar system compared to distribution of captured bodies by Earth-Moon & Jupiter using vo = 8km/s. 
We find that the distribution of captured bodies tends to be more eccentric and highly inclined compared to the distribution of known small bodies in the solar 
system. Long-term evolution of the captured population in necessary before making any concrete predictions. 


4.4 Long term evolution and survival of captured bodies 


We perform long term evolution of a subset of captured objects in 
NEO orbits with a < 20 au. The simulations are carried out until 
a duration of 10 Myr to understand how many objects survive over 
this duration. In total, we perform ~ 20k long-term simulations. The 
orbital parameters of the captured objects is stored for every 1000 
years of evolution. 


We plot the survival fraction as a function of the lifetime of the 
captured interstellar NEOs in Figure 8. We note that the survival 
fraction only accounts for survival in NEO orbits. Thus, objects that 
are still bound but diffused into orbits with larger perihelion and 
semi-major axis are not taken into account. 


The survival fraction decreases monotonically and only 0.1-1% 
of the captured objects survive in NEO orbits by 10 Myr. This is 


significantly shorter than the median lifetime of known NEOs, which 
is around 10 Myr. The captured objects, on account of being highly 
eccentric and prograde, tend to survive for shorter durations. Some 
objects collide with the planets or the sun during this time while other 
objects raise their perihelion upon close encounters with planets and 
move out of NEO orbits. 


Unlike Napier et al. (2021b), we stop our simulations at 10 Myr 
since we are mostly interested in understanding the characteristic 
lifetime of these captured objects. Due to the shorter duration, it is 
hard to ascertain what the functional form of the survival fraction 
is as a function of the lifetime. We do notice that at t > 1 Myr, 
the survival fraction approaches a ~ t~!-© form for both objects 
captured by Earth and Jupiter. This is, in principle, similar to the 
findings of Napier et al. (2021b) although we would need longer 
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Figure 8. Survival fraction of captured interstellar NEOs as a function of 
the lifetime in bound NEO orbits. We notice that the survival fraction mono- 
tonically decreases and only 0.1-1% of the captured objects survive in NEO 
orbits by 10 Myr. This is significantly shorter than the lifetime of known 
NEOs. Once leaving the NEO orbit, the object either becomes unbound or 
switches to an orbit with a larger perihelion and semi-major axis. Objects 
captured by Earth on average survive 2-3x longer than objects captured by 
Jupiter. 


duration simulations to confirm this which are outside the scope of 
the current work. 

We denote the time at which the surviving fraction is half the 
original fraction as T 9.5 and the time where the surviving fraction is 
one-tenth of the original fraction as t9,;. Calculating To.5 from the 
data we find that 79.5, jup ~ 0.05 Myr whereas 10.5 earth ~ 0.13 
Myr. We find 79.1, jup ~ 0.8 Myr and 79.1 earth © 2-1 Myr. This is 
interesting since it turns out objects captured by Earth survive longer. 
The survival fraction for captured ISOs by Earth is almost 10x larger 
than that of objects captured by Jupiter at 10 Myr. We hypothesize that 
this might occur due to a larger fraction of Jupiter captured objects 
on chaotic, often retrograde orbits but a more thorough analysis is 
required which is beyond the scope of this study. 

The objects that are able to raise their perihelion beyond NEO 
orbits survive much longer. We noticed that the overall survival rate 
increases by as much as 10x when all bound objects, regardless 
of whether they survive in NEO orbits, are taken into account. In 
that case, Tt .5 for all captured objects mirrors the results of Napier 
et al. (2021b). This is a noteworthy point since the objects that were 
simulated in Napier et al. (2021b) had much larger semi-major axes. 
This suggests that close encounters in the solar system are often able 
to increase the semi-major axes of captured objects quite efficiently 
even for objects that start in really bound orbits of a < 20 au. 

Although not presented here, we notice that the survival fraction is 
independent of initial vo of the ISOs. This is similar to the findings 
of Napier et al. (2021b) where the authors found that the initial 
hyperbolic excess velocity played no role in determining the lifetime 
of captured objects. 
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Figure 9. Velocity averaged cross-section of capture of ISOs ((a-Voo)) as a 
function of the threshold semi-major axis a;. The shaded region denotes the 
uncertainty arising from the uncertainty in the LSR velocities. The capture 
cross-section is presented for captures by both Earth-Moon and Jupiter. As 
noted in previous studies, we expect the capture rate to be dominated by 
Jupiter, even for NEOs. Jupiter is able to inject objects into NEO orbits at a 
rate that is 104x that of Earth. 


5 DISCUSSION 


Using the capture cross section calculated in the previous sections, 
we now proceed to calculate the rate of capture of objects under 
various velocity distributions and object sizes. 


5.1 Capture in present day solar system 


To calculate the rate of capture, we first calculate the velocity aver- 
aged cross-section given as, 


(OV) = J fo)olwvdmadve (10) 
0 


where f (vo) is the distribution function of the hyperbolic excess ve- 
locity, 0 (væ) is the functional form of the cross-section as a function 
of væ. We assumed that the distribution of excess hyperbolic veloci- 
ties is isotropic. We are interested in understanding rates of captures 
of objects due to close encounters with the Earth-Moon system and 
NEOs arising out of captures by Jupiter. To model the functional 
form for o (væ), we used an interpolation routine rather than assum- 
ing a fitted form. For values of vo > 15 km/s for Earth-Moon and 
Væ > 20 km/s for Jupiter, we assumed © = 0 since we did not obtain 
any datapoints beyond these væ values. Thus, our calculated values 
of capture rates represent the lower bound. 

In present day, we expect the hyperbolic excess velocity distri- 
bution of objects to resemble the field stellar velocity distribution. 
Accordingly, to estimate the rate of capture in present day solar sys- 
tem, we use the known values of the Sun’s velocity relative to the 
Local Standard of Rest (LSR). The solar velocity relative to the LSR 
is given as (vp), vy, Vy) = (1041, 1142, 740.5) km/s (e.g., Schön- 
rich et al. 2010; Tian et al. 2015; Bland-Hawthorn & Gerhard 2016). 
The velocity dispersion of the local stars around the LSR is given 
as (oy, ov, ow) = (33 + 4,38 + 4, 23 + 2) km/s (e.g., Anguiano 
et al. 2018). Using the method described in Siraj & Loeb (2020), 
we construct the probability distribution of incoming speeds of ISOs 
with respect to the Sun. For more information, we refer the interested 
reader to Siraj & Loeb (2020). We note that we assumed the speeds 
are distributed isotropically, an assumption that might not hold in 
reality. 


We plot the velocity averaged cross section ({@Voo)) as a function 
of the maximum semi-major axis of captured objects (a;) in Figure 
9. We find that (ovo) jup/(7Voo)earth ~ 104. Thus the majority of 
the captured interstellar NEOs is due to Jupiter rather than Earth. It is 
a monotonically increasing function of a; indicating that captures at 
larger semi-major axes are easier than captures at smaller semi-major 
axes. Note that this is only for captures due to close encounters with 
the planets. Wide encounters do not play a role in capturing objects 
with a < 100 au for both Earth and Jupiter and can be safely ignored. 

Our results indicate that in the present day, capture is hard due 
to the relatively large velocity dispersion of field stars with respect 
to the Sun. A caveat of our results is the assumption that incoming 
ISO velocities reflect the field star velocities. However, gravitational 
kicks from planets lying in the habitable zone of M dwarf stars (the 
most abundant type of star) can eject objects at ve ~ 50 km/s (Siraj & 
Loeb 2020), where ve is the ejection velocity. If we assume that the 
ejection velocity is 50 km/s, we find that (a7 voo) is almost 4x lower 
than that from before. This indicates the sensitivity of the results to 
the origin of ISOs and highlights the importance of modelling the 
origin of these objects properly. 

Capture would have been possible in the solar birth cluster as well. 
Although we don’t present a thorough analysis of capture in the birth 
cluster, assuming a functional form similar to that used in equation 34 
in Napier et al. (2021a) allows us to derive the approximate rates for 
capture in the solar birth cluster. The functional form for the ejection 
velocity is derived from Moorhead & Adams (2005) and reflects the 
ejection due to a giant planet like Jupiter. Doing so, we find that 
(7Voo)jup (ar = 10au) ~ 10~2au3yr7! which is almost 10°x that 
in the present day environment! However, the configuration of the 
solar system in the birth cluster is a debated topic and the presence of 
nearby stars and the tidal field of the cluster would affect the capture 
and retention process. In light of these uncertainties, we defer a 
careful analysis of capture in the birth cluster to future studies. 


5.2 Present day population and prospects for detection 


To understand the present day population of interstellar NEOs, we 
calculate the rate of capture using the results from the previous 
section. The rate of capture (I) is defined as follows: 


T = nol Vo} (11) 


where nygso is the number density of ISOs. Using 1/‘Oumuamua, 
the number density of Oumuamua sized objects has been determined 
to be 0.1 — 0.2au73 (e.g., Do et al. 2018; Moro-Martin & Norman 
2022). To estimate the size distribution of objects (njgQ(d)), we use 
the following normalization from Siraj & Loeb (2022b): 


d\'4 
niso(d) = nIs0,0 Fal (12) 
(0) 

where njso,o is the number density of ‘Oumuamua sized objects, 
d is the diameter of the object, dg is the diameter of “Oumuamua 
sized objects, and q is a power-law index where 4 < q < 4.5. To be 
conservative, we set njs0,0 = 0.lau~3 do = 100m. This equation 
is valid for object sizes smaller than that of “‘Oumuamua and we find 
that the number density obtained from the normalization used above 
agree with those obtained from other works using double power law 
models (Moro-Martin & Norman 2022). 

The number of ISOs (Niso) in NEO orbits can then be estimated 
as 


Niso = Tt (13) 


where T is a measure of the lifetime of such objects. To simplify 
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our calculations, we assume T = 79.5, 70,1 Obtained in the results 
section. Note that r might vary depending on the semi-major axis 
of the captured objects but our assumption is valid in giving us an 
order-of-magnitude estimate on the current number of captured ISOs 
amongst NEOs.. 

We plot the estimated number of captured interstellar NEOs (Niso) 
of different sizes by Earth and Jupiter as a function of the threshold 
semi-major axis a; using T = 79.5 in Figure 10. Small interstellar 
objects of sizes 1/100 that of ‘Oumumua (~ 1 m) would be captured 
more efficiently owing to the large number density. We note that 
according to our calculations, a steady state population of > 104 ~ 1 
m sized objects should be present with a < 5 au among present day 
NEOs. The number increases to 5 x 10° for a < 10 au. Our results 
also indicate that, although fewer, we would still expect a steady 
state population of ~ 100 0.1do (~ 10 m) sized objects among 
NEOs today with a < 10 au. The steady state population drops 
fast with increasing radius and we find that only 0.2 objects of size 
0.5do (~ 50 m) should exist among NEOs today and only < 0.01 
‘Oumumua sized objects in NEO today. 

The volume capture rate for captured ISOs was calculated to be 
0.051au7yr~! in Hands & Dehnen (2020). Only 0.033% of the total 
captured population was found to be lying within a < 6 au. We find 
that using the same nyjso and T as above, it would imply a Njsg =% 300 
for ~ 10 m sized objects with a < 6 au. This is consistent with the 
Njso for similarly sized objects in our study indicating that our study 
is consistent with previous works. 

Since the number is highly sensitive to the lifetime of these cap- 
tured objects, we perform a similar analysis where we use T = 79,1. 
We present a similar analysis comparing the number of captured 
ISOs in the previous and the optimistic scenarios in Figure 11. This 
represents a more optimistic scenario. We find that the Niso in this 
case is approximately 10x larger than that of before. In this case, we 
expect to find a steady state population of 2-5 ~ 50 m sized objects 
with a < 10 au among NEOs today . This has major implications for 
the future detectability of such objects as we explain below. 

We consider the detectability of such objects using the Vera Rubin 
Observatory (LSST). LSST has a limiting magnitude of ~ 24.5 ina30 
second exposure. Assuming an albedo similar to that of ‘Oumuamua, 
we find that the absolute magnitude of 50 — 60 m sized objects in 
NEO should be around the limiting magnitude of LSST. Assuming 
these objects have a ~ 10 au would imply a period of rotation of 
~ 30 years. This would imply that over a 10 year run of LSST, we 
should be able to detect ~ Nīso Ta = L Niso. Using an optimistic 
estimate for Nīso implies that LSST should be able to detect O (1) 
such captured interstellar NEOs during its lifetime. If these captured 
objects had larger albedo, LSST would be able to detect even smaller 
objects going down to 20 — 30 m sized objects enabling O(10) 
detections over its lifetime. Similar analysis can also be performed 
with the Hyper Suprime Cam (HSC) which gives us prospects of 
detectability close to those obtained using LSST. 

It is evident from our study that smaller captured ISOs outnum- 
ber larger ISOs. However, detecting these 1 — 10 m sized objects 
is a hard task. Shao et al. (2014) propose synthetic tracking to de- 
tect very small objects down to sizes of 7 m among NEOs. Such a 
method would allow detections of smaller interstellar NEOs and en- 
able O(10) detections per year. Machine-learning (ML) based meth- 
ods (e.g., Hefele et al. 2020) may also improve the ability to detect 
smaller objects. We plan on investigating ML based techniques for 
detection in the future. Prospects for detectability using future space 
based observatories such as NEO Surveyor should also be investi- 
gated in the future. The detections should focus on objects with orbits 
that are largely prograde and highly eccentric although work needs 
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to be done to properly ascertain the distribution of orbital parameters 
after the object has undergone encounters with planets in the solar 
system after getting captured. This is beyond the current scope of the 
study and will be investigated in the future. 


6 CONCLUSIONS 


Interstellar objects present a unique mechanism to investigate the 
formation and evolution of planetary systems including our own. 
Although rare, ISOs can be captured into bound orbits by different 
planets in the solar system. Whether any captured ISOs exist in th 
present-day solar system is a major area of interest. A population of 
captured ISOs can provide valuable information regarding the origins 
of such objects. 

In this study we examined the capture of ISOs into near Earth orbits 
using a large suite of N-body scattering simulations using a new 
hybrid integration scheme. We computed the capture cross-section 
of Earth-Moon and Jupiter and found that Jupiter should dominate 
the capture of interstellar objects into near Earth orbits by a factor of 
104 compared to that of Earth-Moon. These captured objects have 
larger semi-major axes, lower inclinations, and higher eccentricities 
than the current known population of Near Earth Objects. We find 
that most captured interstellar NEOs in the solar system arise from 
beyond the orbit of Saturn. 

Long term simulations of such objects reveals that they’re more 
chaotic than known NEOs and have shorter lifespan in NEO orbits 
than the known population. The half-life of the Jupiter captured pop- 
ulation is ~ 0.05 Myr whereas that of the Earth captured population 
is x 0.13 Myr. It takes ~ 0.8 Myr for the Jupiter captured population 
to dwindle to 10% of the original fraction whereas it takes 2.1 Myr 
for the Earth-captured population to do the same. 

The computed cross section is then used to determine the capture 
rate in the present day solar system assuming the ISOs originate with 
a velocity distribution similar to that found for the field stars in the 
solar neighborhood. The capture rate is computed for objects of dif- 
ferent sizes. Using the computed life time of survival, an approximate 
number of ISOs among NEOs in the steady state is determined. 

We find that smaller (~ 1 m) sized objects should dominate the 
captured population and in the optimistic scenario, can make up a 
population of 10° objects with a < 10 in NEO orbit. ~ 10 m sized 
objects should be present in the range of 10-1000 and ~ 50 m sized 
objects should be present in range of 0-6 with similar orbital param- 
eters. Unfortunately, according to our simulations, no Oumuamua 
sized objects are able to be captured and retained in NEO orbit. 
However, such objects may exist within the solar system but not in 
NEO orbits. 

The captured population can have major implications for de- 
tectability. Although LSST would only be able to resolve ~ 50 m 
sized objects with a < 10 au, other mechanisms such as synthetic 
tracking should be able to detect objects as small as 7 m among 
NEOs. Surveys should focus on objects with orbits that are highly 
prograde, eccentric, and have semi-major axes of a ~ 10 au. How- 
ever, to better ascertain the distribution of orbital parameters, long 
term simulations should be used. 

To our knowledge, this is the first study that explores the possibility 
of presence of interstellar objects in our neighborhood among NEOs. 
Interstellar NEOs would be preferable to other captured ISOs due to 
the proximity to Earth which would aid in detectability and possible 
visits by space probes. With a substantial number of potentially de- 
tectable objects hiding among NEOs, we demonstrate the need and 
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merit for more research in the areas of dynamics and detection of 
these small-sized objects in the future. 
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Figure 10. Number of ISOs (Niso) present in steady state as a function of the threshold semi-major axis a;. The shaded regions represent the numbers calculated 
using different power law index parameters as given in equation 12. The lower limit is taken as q = 4 and the upper limit is taken as q = 4.5. Left: captured 
ISOs as a result of close encounters with Earth-Moon. Right: captured ISOs as a result of close encounters with Jupiter. Smaller sized objects (~ 1 m) should 
outnumber larger sized objects. We expect O(10°) of these objects to lie within the orbit of Saturn and in NEO orbits. ~ 10 m sized objects number are also 
present in NEO orbits today, according to our calculations. We expect to find a steady state population of cov” such objects. This calculation represents the 
conservative estimate calculated using a survival lifetime of T = To,5. Under this conservative estimate, we should not expect to find any Oumumua sized objects 


among NEOs today. 


Period [yr] 
89 354 1000 


32 


101 102 
ar [au] 


Figure 11. Comparison of the number of captured ISOs in NEO assuming the original and optimistic scenario. The optimistic scenario is denoted in the darker 
shade. The colors denote the size of the object and are the same as that in Figure 10. An optimistic scenario predicts almost 10x as many captured ISOs compared 
to the original scenario. In this case, we predict a steady state population of 3-6 ~ 50 m in NEO with a < 10 au. Ina 10 year lifespan of LSST, we can assume 
to detect ~ 1 such object assuming that the brightness profile is similar to that of Oumuamua. 
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APPENDIX A: COMPARISON TO IAS15 


In order to verify the validity of our hybrid integration scheme, we run 
a set of ~ 10’ scattering experiments with the Earth-Moon system 
using the IAS15 integrator (Rein & Spiegel 2015) and compare the 
orbital parameters of the captured objects to those obtained using our 
hybrid scheme. For the sake of clarity, we only present the results 
from væ = 8 km/s. We plot the distribution of semi-major axis (a) 
and inclination (i) obtained from using IAS15 and our hybrid method 
in Figure Al. We notice that the distributions of captured orbital 
parameters is statistically similar using both methods. The median 
semi-major axis of the captured objects using IAS 15 is 44.6 au while 
that using our hybrid method is 46.2 au. The 10 and 90 percentile 
values of the captured semi-major axis using IAS15 are 13.1 au and 
299.8 au respectively. Using our hybrid method, we obtain 13.2 au 
and 322.8 au for the 10 and 90 percentile values of the captured 
semi-major axis respectively. The median inclination using I[AS15 
was found to be 11° whereas that using our hybrid method was found 
to be 10.5°. This indicates that the results obtained using our method 
are statistically similar to those using IAS 15, confirming the validity 
of our simulations. Although not presented here, we performed a 
similar analysis with smaller values of væ and found similar results 
thereby confirming that our method is able to simulate both weak 
and close encounters properly. 
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Figure A1. Distribution of the orbital parameters of captured objects for vo = 
8 km/s using the IAS15 integrator and our hybrid method. We notice that the 
distributions of semi-major axis (a) and inclination (i) of the captured objects 
are very similar between the two schemes. This validates the robustness of 
our hybrid method. 


